Persistent currents of noninteracting electrons 
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We thoroughly study the persistent current of noninteracting electrons in one, two, and three 
dimensional thin rings. We find that the results for noninteracting electrons are more relevant for 
individual mesoscopic rings than hitherto appreciated. The current is averaged over all configura- 
tions of the disorder, whose amount is varied from zero up to the diffusive limit, keeping the product 
of the Fermi wave number and the ring's circumference constant. Results are given as functions of 
disorder and aspect ratios of the ring. The magnitude of the disorder-averaged current may be larger 
than the root-mean-square fluctuations of the current from sample to sample even when the mean 
free path is smaller, but not too small, than the circumference of the ring. Then a measurement 
of the persistent current of a typical sample will be dominated by the magnitude of the disorder 
averaged current. 

PACS numbers: 73.23.Ra, 73.21.-b 



I. INTRODUCTION 

One of the consequences of the Aharonov-Bohm (AB) 
effeclji is that a finite normal (i.e. non-superconducting) 
mesoscopic ring exhibits a persistent current (PC) when 
the AB magnetic flux through its opening is non zeroi^^— 
The PC does not decay with time when the dephasing 
and the thermal lengths are larger than the ring cir- 
cumference. This results from the fact that the PC re- 
flects an equilibrium state even when the ring has a fi- 
nite resistance due to defect scatteringj^i^ii The PC is 
periodic in the flux $ with a period given by the mag- 
netic flux quantum $o = 2TThc/e. Measurements of the 
PC^-ii often stimulated the theoretical studiesJ^".2i To- 
day, this fundamental phenomenon of quantum mechan- 
ics still challenges both theoreticians and experimental- 
ists of mesoscopic physicsi^"— Persistent currents are 
also relevant for the orbital response of semimetals and 
aromatic molecules^ and for the ongoing interest in 
nanotubesj^ 

At zero disorder, the azimuthal component of the ve- 
locity associated with each single-particle eigenstate of 
the Hamiltonian of noninteracting particles is shifted due 
to the AB flux $ < $o/2, by Ai; = 2nh(j)/ML. Here M 
is the electron mass, L is the circumference of the ring 
and (/)= $/$o- One may naively assume that the current 
density is —neAv, where n is the density of the electrons. 
In a normal ring, because of level crossing, the occupa- 
tion of the levels changes with the flux. As a result, once 
level-crossing occurs, the PC density of the normal ring 
is much smaller than —neAv. In a superconducting ring 
—neAv gives the value of the PC density at zero tem- 
perature and zero disorder. It might be argued that in 
a perfect superconductor at zero temperature, the above 
occupation switching is suppressed. Thus, the attractive 
interaction in a superconductor, which enforces the pair- 
ing correlations, strongly enhances the PC compared to 



the normal state value. Note that the current of a su- 
perconducting ring is an intensive quantity -it does not 
depend on the size of the system. In the normal state, 
the current is only a mesoscopic effect-proportional to an 
inverse power (-1 in the ballistic ID case^) of the system's 
length. 

The current of noninteracting electrons in 2D cylinders 
in the grand-canonical ensemble was studied analytically 
in the limit of zero disorder and in the diffusive limit J^"— 
In these works the PC was calculated in two geometries: 
"short" cylinders, H <^ L, where H is the height of the 
cylinder, and "long" cylinders, H ^ L. Cheung et alt^ 
studied the case of a 3D short and thin diffusive cylinder 
as well. In the zero-disorder limit, the PC was calcu- 
lated by summing the velocities, with appropriate fac- 
tors, of all the states that, after the energy shift due to 
the flux, are below the Fermi energyJ^ In the diffusive 
limit the PC may be averaged over the configurations 
of the impurities. It can be calculated as a function of 
the magnetic flux from the density of states in the dif- 
fusive limit li^iii Entin-Wohlman and Gefeoi^ calculated 
the impurity-ensemble-averaged current of long cylinders 
using the linear response theory in (p, which is valid only 
for (/)< 1/2. 

Our work extends the above researchi^^— in two ways. 
First, we describe the current for any degree of disorder 
between the previously studied limits of perfectly clean 
systems and diffusive systems. Second, we consider 3D 
thin rings with a flnite width W for which W L {in 
contrast to IF < a, where a is the smallest microscopic 
length of the system) F^^iiiii^ We also correct, and gen- 
eralize for any given value of the flux, the expression for 
the PC as calculated by Entin-Wohlman and Gefen for 
"long" 2D cylinderSfi^ In the latter, a calculation error- 
gave a result of incorrect sign and magnitude for the pref- 
actor of the dominant (for L ^ £, where £ is the elastic 
mean free path) exponential dependence. 
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The expressioni^ for the disorder-averaged PC in the 
grand-canonical ensemble at zero temperature is given in 
SecHH This expression can be simplified in two regimes, 
defined in Sec. Illlj which we name the uncorrclated and 
the correlated-channel regimesi^ In sections ITVl and IVl we 
preform the simplifying steps that are allowed in each 
regime, and then obtain the leading-order expressions 
of the PC in the zero-disorder and the diffusive limits. 
The specific conditions for which these two limits hold in 
both the uncorrelated and the correlated-channel regimes 
are given in Table HI In Refs. [isl - fTsI the same simplify- 
ing assumptions had been used, but were referred to as 
"short" and "long" cylinders. We find that these picto- 
rial definitions do not agree with the regimes in which 
the corresponding results hold. Our results for PC of 2D 
cylinders in the zero-disorder and the diffusive limits for 
the uncorrelated-channel regime, and in the zero-disorder 
limit for the correlated-channel regime, agree with the 
ones obtained in Refs. [Tl-[ll For 2D cylinders, our re- 
sult for the PC in the correlated-channel regime in the 
diffusive limit is new. 

The disorder-averaged PC is highly sensitive to the ex- 
act value of kpL, as it contains a factor of sin(fci?i), 
where kp is the Fermi wave number. In Sec. IVII we dis- 
cuss the way to compare the measured average PC in an 
ensemble of rings to the theoretical results depending on 
the variance of the value of kpL among the rings. In 
this section the disorder-averaged PC is also compared 
with the root-mean-square (rms) fluctuation o^^'^" of the 
PC with respect to the disorder. We find that as long 
as the system is not too diffusive, the magnitude of the 
disorder-averaged current may be larger than the current 
rms fluctuations. As discussed in Sec. IVID our result 
for the disorder-averaged PC of noninteracting electrons 
agrees with the PC measured in a 2D clean annulus by 
Mailly et al^, but has a larger magnitude than the one 
measured by Rabaud et alJ^ The results of our study are 
discussed in Sec. I Villi 

In contrast with the Green function technique used in 
the main body of this paper, we give in the Appendix □ 
a novel approximation for the PC of a 3D ring in the 
zero-disorder limit. This approximation is based on the 
canonical ensemble results for a ID ring, and on the prob- 
abilities that, at a given flux, the number of electrons in 
a given transverse channel is odd or even. 



The AB flux, which does not penetrate the ring itself, 
is given by the magnetic vector potential A = (^$/27rr, 
where r is the radial coordinate and is a unit vector ori- 
ented along the ring. The eigenstates of "H, in cylindrical 
coordinates, are 



^(r,^,z)=e"'^sin(^) 



where n — 0, ±1, ±2, 



I, 2, .., and 



k = ^2Me/h? - [Tiq/HY . 



(2) 



(3) 



Here J and Y are the Besscl functions of the first and 
second kind. The boundary conditions 'ip[r = (L/27r — 
W/2)] = i}j[r = (i/27r + W/2)] = set the ratio be- 
tween the prefactors Ci and C2 and the eigenenergies. 
For W the eigenenergies are given hy^ 



f 



1 S^-l [2(71 + 0)] 



L2 



O [{W/Lf 



(4) 



where s is a positive integer. In this work, all energies 
arc shifted so that the single-particle ground state energy, 
for which q = s = l^n ~ (p ~ Q, is zero. We henceforth 
neglect the term of order {W/L)'^ in Eq. (g]). 

We now introduce disorder, induced by impurities 
having point-like potentials. The PC, averaged over a 
grand-canonical ensemble of disordered systems having 
the same mean free path but different impurity configu- 
rations, is given byi^ 



X [G+{[q,s,n],E)^G-{[q,sM,E)\ 4°' . (5) 

Here the Fermi distribution function, f{E), sets the 
chemical potential as an upper bound on the integra- 
tion at zero temperature. The current associated with a 
single-electron wave function is given by 



2'iThe , 



(6) 



II. THE EXPRESSION FOR THE PERSISTENT 
CURRENT 

In this section we obtain an expressioni^ for the 
impurity-ensemble-average zero-temperature PC of non- 
interacting electrons. We consider spinless electrons in a 
ring of a mean circumference L, a width W ^ and a height 
H. In the absence of disorder, the Hamiltonian is given 
by 



(1) 



In Eq. ([S]) , the disorder-averaged retarded and advanced 
Green functions are denoted by and G~, respectively. 
The expressions for the disorder-averaged Green func- 
tions, for kpl ^ I and within the Born approximation. 



G^{[q,sM.E) 



E 



ih 



(7) 



where t is the elastic mean free time. Equation ([5]) for 
the disorder-averaged PC is given as a sum over channels 
(q,s). However, in the corresponding expression for the 
non-averaged current, one should use the non-averaged 
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Green functions and consequently for a specific configu- 
ration, the channels are mixed in the expression for the 

We note that the (q, s) term in Eq. ([5]) is given by 
the averaged PC in a ID ringii with a shifted chemical 
potential 

fi fi{q, s) = e{q, s, n = 0, = 0) , (8) 

namely, 

{I)=Y^{I^^[^,{q,s)]) . (9) 
The current of a ID ring, calculated in Ref. [l3, is 

(/i^)=2/o> — ^ cos ( mk F L) e" ^ . (10 

^ ' 7rm 

m—l 



Here Iq = evp/L , where uj;^ is the Fermi velocity)^ 
In Eq. ([9]) each (q, s) term has its Fermi wavenumber 
determined by Eq. dS]) 



kF{q,s) = kF^/Ti<hs)]Ti ■ (11) 

Equation (ITOl) is valid for /i ^ {/i/2r, a}, where a = 
l-K^f? jMl? is the prefactor of (n-|-(/))^ in the expression 
for the eigenenergies, see Eq. (j4]). 

Substituting the ID result, Eq. ([13 , in Eq. (O, we 
obtain that at zero temperature 

oo 

(/) = (^'») sin(2^m0) , (12) 

m— 1 

where the disorder-averaged harmonics are given by 



g— 1 s— 1 ^ v-i // / 



The approximate numbers of the occupied channels cor- 
responding to momenta in the radial and the z directions 
are 

= fc^VF/TT , iV^ = kFHj-K , (14) 

respectively. In the upper bounds on the summations 
over q and s, one needs to take the closest integer values 
for iVf. and from below (but not less than one). 

In Eq. ([13]) we sum over the contributions of the occu- 
pied channels, which obey {s / + {q/ N < 1, so that 
li{q, s) > 0. In a diffusive system, one might worry about 
the contribution to (Im) of channels with high transverse 
momentum which satisfy 

£[kF{q,s)/kF]<l/kF{q,s) , (15) 

and are therefore not diffusive. Their contribution is 
given by an expression similar to Eq. (jlOp . where a term 
of s/AkpL multiplies the exponent and divides Iq. In 
Eq. (jl3p we ignore this extra reduction, since only a few 
channels may satisfy Eq. (|15p and their contribution to 
the PC is anyhow small. 

III. APPROXIMATIONS FOR THE PC 
HARMONICS 

In this section we identify different regimes in which 
the expression for the disorder-averaged harmonics, see 
Eq. can be simplified. 



A. Dimensionality of the system 

The ring is considered to have a significant thickness 
along the radial direction when Nr ^ I [see Eq. (|14p ] 
and when the ratio between the exponential in Eq. (|13p 
with a small index s to the following s + 1 term is much 
smaller than, say, 10. Thus, for the calculation of (Im) 
many s values give significant contributions when 

kpW , P^] 

^»|land^-|. (16) 

When the "much larger" sign in Eq. (|16p is replaced by a 
"smaller" or "comparable" one, the ring is considered to 
be of zero dimension along the radial direction, and we 
use only s = 1. 

Note that condition ([T5)) depends on L/£. This can be 
understood by the following argument: The phase of the 
Green function of a particle that encircles the ring de- 
pends on the specific disorder configuration. Averaging 
the PC over all configurations of disorder results in the 
exponential decay of (/^^), see Eq. ((TU)) >^ In a multi- 
channel ring, the overall path, and correspondingly the 
variance of the phase shifts, increase as the transverse 
momentum increases. This results in the increase of the 
exponential decay rate in Eq. (jl3p for increasing chan- 
nel index. Indeed, as we see in Eq. (fT6)) . increasing the 
disorder may decrease the effective dimensionality of the 
system. The condition for considering the ring to have 
a finite height is given by Eq. (jl6p upon replacing W 
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with H. In this way the system is classified as one of 
the following: ID, 2D annulus, 2D hollow cylinder, or a 
3D ring. In the 2D annulus case one sums over s taking 
q — 1^ and in the 2D cylinder case the summation is over 
q keeping s = 1. 



we use the approximation 



cos [mkp{q, s)L] cos [mfci^(g', s)L] = , 



and obtain 



(18) 



B. Contributions of consecutive channels to (Im) 

The discrete summation over the channel indices in 
Eq. (|T3|) makes the expression for (/„) hard to handle an- 
alytically. In this subsection we define two regimes where 
one can overcome this difficulty. The contributions to the 
mth harmonic of consecutive transverse channels (s and 
s + 1, or g and q + 1) are uncorrelated when the change 
in the arguments of the corresponding cosine terms, see 
Eq. (fT5|) . is larger than, say, 7r/4. This difference between 
the arguments of the cosines increases with increasing 
channel index. Hence, if the lowest two transverse in- 
dices obey this condition, then higher indices will fulfill 
it as well, so that all channels arc uncorrelated. Thus, the 
channels associated with the z direction are uncorrelated 
when 



H 
mL 



< 



2-K 
kpU 



(17) 



( {Imf ) 



1/2 y/2 



kp 



cxp 



mL 



ikpiq, s)/kf 



(19) 



We have confirmed numerically that the standard devia- 
tion of {Im) obtained from Eq. ((T3)) gives the same value 

for ( (/m)^ )^^^ as given by Eq. ([TO]) . For the calcula- 
tion of the standard deviation of {Im) wc have inserted 
in Eq. (|13p the parameters of the ring used by Mailly et 
al.^ see Sec. IVIIi and considered many values of kpL in 
a segment of a width of IOtt. 

When the first harmonic is in the uncorrelated-channel 
regime, the harmonics with m up to m ^ SkpW'^i/Tr'^L 
are also in that regime, see Eq. ((T6)) . In this case, the 
contribution of higher harmonics is negligible. Therefore, 
in the approximate expression 



The same rule applies to channels of consecutive s in- 
dices upon replacing H with W. The regime defined by 
Eqs. ((T5)) and (|17p will be referred to as the uncorrelated- 
channel regime. 

In the uncorrelated-channel regime the dependence of 
the PC on the parameters kpL, N^, and iV^, which appear 
in the arguments of the cosines in Eq. (jl3p . is non-trivial. 
This is demonstrated in Fig.[TJ We thus turn to calculate 
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FIG. 1: The disorder-averaged PC depends on kFL,Nz and 
Nr in a non-trivial fashion. Here we plot (l^^i), seeEq. 
for L/^ = 5 and for iV^ = 70 (solid line) and N^, = 100 (dashed 
line). The typical magnitude of the disorder- averaged current, 
given by Eq. (|19|) . for the above two values of is 0.25/o 
and O.SO/o, respectively. 

the typical magnitude of the disorder averaged harmonics 

( {Im)'^ Yl"^ . The overline denotes averaging over kpL 
within a segment 8(kpV) <^ kpL of a width of > 27r. 
Note the different notations of averaging over kpL and 

averaging over disorder. In the calculation of ( (/„i)^ )^^^ 



1/2 



\ m=l 



sin2(27rm0) , (20) 



we can use the expression given in Eq. (|19p for {Im)" for 
all the relevant harmonics. 

For a 2D cylinder, the maximal q whose contribution 
to (/„j) is not negligible, see Eq. (|13p . is 



(21) 



When Eq. p6)) is satisfied and the cosines of sequential 
indices with q < (j™a,x are correlated, then the sum in 
Eq. (fT3l) can be replaced by an integral. Since the differ- 
ence between the arguments associated with sequential 
channels increases as the index of the channel increases, 
the condition for the channels to be correlated is 



mL [fc^^(C-l,l)-fc;^(Cax,l)] < 



(22) 



N, , the condition 



has the form 



When g™^^ = 
H/L > lOm^kpL. The correlated-channcl regime for a 
2D annulus is defined in the same way, but the limitation 
W L oi our analysis makes this regime irrelevant for 
that geometry. We refer to this point in more detail at 
the end of Sec. El The expressions for the conditions for 
the uncorrelated and the correlated-channel regimes, in 
the zero-disorder and the diffusive limits are summarized 
in Table m 
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Conditions associated with the z direction 


Results 


uncorrelated: zero disorder 
uncorrelated: difTusive 
correlated: zero disorder 
correlated: difTusive 


1 < kpHI-K < 2L/H 
^L/8£ < kpH/iv < 2L/H 
H/L > lOm^kpL 
H/L > max{WO£/ H,TT/kFL} 


(|24I25P [3D rinesj 
(I2S1) [3D rings] 
(f29)) [2D cylinder] 
(|30l) [2D cylinder] 



TABLE I: The results for the PC in the zero-disorder {£/L — >■ cxd) and the diffusive {£ <^ L) limits. The conditions defining 
the uncorrelated and the correlated-channel regimes are given in the second column for a 2D cylinder. For a 2D annulus the 
conditions are the same with H replaced by W . For a 3D ring the conditions should be satisfied for both azimuthal directions. 
In the third column we refer to the appropriate expressions for the PC. 



IV. UNCORRELATED-CHANNEL REGIME 

Consider a 3D ring in the uncorrelated-channel regime, 

defined by Eqs. (US]) and To estimate ( {Im)'^ f/"^ 

we replace the sum in Eq. p9p by an integral over x = 

Ntot = 



j)2 + (s/Nr)^, and add the factor 2a;A'tot, where 
jNrNz is the total number of occupied channels 



1/2 2^2 




(23) 



In Fig. [2] the magnitudes of the first and second harmon- 
ics are plotted as a function oi L/£ using Eq. (P^ . Here 
one can see that with increasing disorder, the first har- 
monic becomes more dominant. 




FIG. 2: The PC of a 3D ring in the uncorrelated-channel 
regime. The typical magnitudes of the first harmonic (solid 
line) and the second harmonic (dash-dotted line) are plot- 
ted in units of loy/Ntot, using Eq. (|23[) . In the inset 

)^^^/Io (solid line) is obtained by substituting 
(|23|) . For a later discussion, the rms flue- 



N^^ Nr = 20 in Eq. 
tuations of the PC with respect to the disorder, 51 /Iq, (dashed 

line) are plotted using Eq. (|33]). Here ( {Im^if = SI 
when L/£ = 8.5, in agrement with Eq. (|34p . The horizontal 
axis of the inset begins at L/£ — 3 since Eq. (|33|l is valid only 
in the diffusive limit. 



From Eqs. jM]) and we obtain^ 

((/3^)')'^'=/ov^Vl0l(l-2|</.|). (25) 

Note the enhancement of the PC magnitude by the 
square root of the channel number. Deep enough in the 
diffusive limit, L/£ > 10 , the PC is dominated only by 
its first harmonic. Here, the magnitude of the PC is given 
by the limit L/^ > 1 of Eq. ^ 



( )'^' = ^y|/oVA^e-^sin(27r. 



0) • (26) 



This reproduces the resulli^ of Ref. [H. 

The PC harmonics of a 2D annulus are given by 



I (l-a;2)exp(^- 



mL 



£Vl - 



dx 



(27) 



Results for a 2D annulus in the uncorrelated-channel 
regime and the zero-disorder limit are given by Eqs. (|24p 
and ([25)1 with A'tot replaced by ANr/S. Here, replac- 
ing Nr with Nz gives the expression for the PC in a 2D 
cylinder obtained^ by Cheung et al— In the diffusive 
limit, the PC of a 2D annulus or a 2D cylinder in the 
uncorrelated-channel regime amounts to mul tiplying the 
expression in Eq. ((26)) by the factor a/ ttL / 8£ and replac- 
ing A^tot by Nr or Nz, respectively. The latter yields 
the results obtained in Refs. [l^ and The difference 
between the powers oi L/£ between the 2D and the 3D 
expressions is due to the difference of the densities of 
states of the transverse channels in these cases. 

The similarity between the PC of a 2D annulus and the 
PC of a 2D cylinder is hardly surprising since these two 
cases of finite width and of finite height are topologically 
equivalent for the AB fiux, and the eigenenergies are the 
same as long a.s W <Si L. 



Equation (j23p can be further approximated in the zero- 
disorder and in the diffusive limits. In the first limit 



V / nm 



(24) 



V. CORRELATED-CHANNEL REGIME 

For a 2D cylinder, the correlated-channel regime is de- 
fined by Eq. ([H]) (with H replacing W) and Eq. ([H]). In 
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this case we replace the summation over q in Eq. p3p by 
an integration and obtain 



cos [mkpL\/ 1 — exp f — 



mL 



2^^/^ 



. (28) 



In the zero-disorder Umit, Eq. ([28|) yields the resulij^ of 
Ref.lil 



j2D _ 



1 



cos(mfcFi- 7r/4) . (29) 



The diffusive limit of the PC of a 2D cylinder in the 
correlated-channel regime is found here to be given by 

, v2 sin(27r(^) l , , , , ^ 

= ^-^=i=^/oiV, e-^ cosikpL - ^/4) . (30) 

(The higher harmonics are negligible.) The conditions for 
the correlated-channel regime in the zero-disorder limit, 
see Table HI cannot be satisfied for the radial direction 
together with the restriction ^ L, for most reason- 
able values of kpL. The limit of a diffusive annulus, see 
Table H is satisfied, for < L, only when L/£ > 130, 
but then the disorder-averaged PC is irrelevant. 

In Fig.|3]the magnitude of the disorder-averaged PC is 
plotted using Eq. pop 0jS &- function of L/t in the dif- 
fusive regime. The results (|29p and ([5(7)) are reduced 



by 1/ \JkpL compared to the results in the uncorrelated- 
channel regime in the zero-disorder and the diffusive lim- 
its, see Sec. IIVI However, these results are enhanced by 
^JWl and by y/Wl{L / 1)^ / ^ , respectively. 
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FIG. 3: The disorder-averaged PC of a 2D diffusive cylinder 
in the correlated-channel regime (solid line) is plotted as a 
function oil/L. We replace sm{kFL) by l/\/2 in Eq. (f30|) to 
obtain the typical magnitude, and use iV^ = 10^ and kpL = 
5x 10^ The rms fluctuations (dashed line), see Eqs. (|31|l and 



(|33|l . equals ^/^j), for the above parameters, at L/£ ~ 10, 
in agreement with Eq. (|35p . Both ^/^^i) and 51 are given in 
units of Jq. 



VI. THE RMS FLUCTUATIONS VERSUS (/) 

The disorder-averaged PC is very sensitive to the exact 
value of kpL, see e.g., the cosine factor in Eq. ([2^ . In 



contrast, the rms fiuctuations of the current in respect to 
the disorde r^^i^" 



51 =[{!-)- {If Y'^ 



(31) 



are not sensitive to kpL. The common practice in PC 
measurements is to determine the total current, /tot, 
from the measurement of the overall magnetic response 
of N rings. This current is related to both the disorder- 
averaged current and to the current rms fiuctuations by 



N (/) ± \/N5I 



±1 



1/2 



±51 



5{kpL) < TT 
5{kpL) > TT 



(32) 



Here 5{kpL) is the variation oi kpL in an ensemble of N 
rings. Equations (|32p hold also for the harmonics (replac- 
ing / by /„i). If the ring is in the uncorrelated-channel 

regime, one may replace (/) by ±( {if )^/^ in the top 
equality of Eqs. ([5^ . while if the ring is in the correlated- 
channel regime, one needs to replace the cosine factor in 

Eq. (EH) for (/,„) by 1/V2 in order to obtain ( {I,nf 
in the bottom equality. 

The rms fluctuation due to the disorder of the h/e 
harmonic of the current for a thin-walled {L ^ {W, H}) 
ring in the diffusive limit is given byi^i^ 



51 



/o sin(27r0) [i < L] 



(33) 



This result is independent of the number of channels, i.e., 
of W and H . These current rms fiuctuations do not ex- 
ist for £/L ^ 1, see Eq. ([51]) . Thus, the contribution to 
/tot which is not related to interactions, is expected to 
be given by Eq. in the zero-disorder limit. Equa- 
tion ([55]) for 51 is strictly valid in the diffusive regime, 
but is expected to give a correct order of magnitude for 
systems in which £ and L are comparable. 

In Figs. [2] and [31 the crossover from the dominance of 
the disorder-averaged PC to the dominance of 51 can be 
observed. In the uncorrelated-channel regime, the typical 
magnitude of the disorder-averaged current of a 3D ring 
is equal to 51 at L/t = 5, 10, 14 for iVtot = 20, 10^, 10^ 
respectively. These values are obtained, for L/i > 1, by 
comparing Eq. with Eq. 1^ 



1/2 



> 61 



(34) 



The analogous result for a 2D cylinder in the correlated- 
channel regime is 



Nz > 0.9- e 

1j 



^'^'^ ^ {{llZ.f)"'>5I. (35) 



For kpL = H/L = 100, the equality ( {l^^^^) f'^ = 51 
is satisfied, see Eq. ([35|), for TVtot = 22, 135, 700 at L/l = 
5, 10, 14, respectively. 
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VII. DISCUSSION OF EXPERIMENTAL DATA 

Since the first harmonic is not expected to be affected 
by electron-electron interactions ; "'^'^i we may compare 
its measurements^"— >2ii^ with calculations of the typical 
magnitude of (/) and SI. 

Mailly et ai— studied the PC in an almoust ballistic an- 
nulus of GaAlAs/GaAs, characterized hy L = 8.5/im, i = 
ll^m,kp = 1.5 X 10^TO~^,u_F = 2.6 x lO'^m/s, and 
W = 0.16/iTO. These parameters, which yield /q = 5nA 
and Nj. ^ 8, satisfy conditions (|16p and ([T7)) for the 
uncorrelated-channel regime. We insert these param- 
eters in our result Eq. (P7)) and in Eq. adding 
a factor of two due to spin degeneracy. This yields 

( (^m=i)^ = 1-44, and SI = 1.3/osin(27r<?!)). We see 

that SI and ( (if^j^i)'^ )^/^ are comparable, and both are 
in fair agreement with the measured PC of (0.8 ± 0.4)/o. 
Using the expression for the PC of a 2D cylinder in the 
zero-disorder limit obtained in Ref. [l^ (replacing H with 
W) yields a value larger by a factor of ~ 2 compared to 
our result. When £ ^ L, the ballistic, diffusive and ex- 
act expressions should give the same order of magnitude 
for the PC. Indeed, using the expression for the PC of a 
diffusive annulus in the uncorrelated-channel regim o^^i^^ 
gives a value that is very close to the one obtained from 
Eq. ([23|) for the parameters of the annulus measured in 
Ref. la 

Rabaud et ali^ measured the PC of an array of 16 
ballistic rings of GaAlAs/GaAs. Those rings are in fact 
squares whose external total edge length is 16/im and the 
internal one is 8fim, yielding L = 12^m. The rings are 
also characterized hy £ = 8iim,kF = 2 x 10^m~^ ,W ~ 
0.8/xm, and vp = 3.2 x lO^m/s, implying Iq = A.2nA and 
Nr = 50. The measured total PC obtained for discon- 
nected rings, divided by the square root of the number of 
ringS)^ was (0.33 ± 0.07)/o. Neither the uncorrelated- 
channel regime nor the correlated-channel regime can 
be associated with these rings, since both Eq. ([T7|) and 
Eq. are not obeyed by the above parameters. There- 
fore, we use our result Eq. (|13p. with q = 1 and a factor 
of two due to spin degeneracy, and obtain values for 
in the regime (— 3/o,3/o), whose standard deviation is 

( {ll?=if Y'^ = l-l^- From Eq. ^ we find that 
SI = 0.7/o sin(27r(/)). The discrepancy between the mea- 
sured value, the above ( (I'^i)'^ )^/^, and SI may be due 
to the geometry (squares instead of rings) as well as due 
to decoherenceFi^ The relative large W may also play a 
role. 

One may compare our result for ( {lj,{^i)'^ )^^^ for 
the parameters of Ref. [13 with results of previous the- 
oretical studies for these "short" annuli.— "— The latter 
correspond to ( (l^^i) )^/^ = 7.5/o in the zero-disorder 

limit, and ( {ify^i)^ Y^"^ = 4.5/o in the diffusive limit 
(as given by Eqs. ((24|) and ([26| . adapted to 2D and in- 
cluding a factor of two due to the spin degree of freedom. 



see Sec. lIVp. Hence, our result is in a smaller disagree- 
ment, compared to results of former studies^"— with 
the measured one. This is due to the fact shown above 
that the conditions for Eqs. ([M)) and (^5)) to be valid are 
not satisfied by the parameters of the rings measured in 
Ref. 0. 

The first harmonic, measured for the diffusive rings 
used in the studies of Jariwala et al.^ and of Bluhm 
et al-r^ fairly agrees with the theoretical value for SI. 
Here the rings are deep enough in the diffusive regime, 

and so ( (/)^ )^/^ ^ SI. In the very recent work of 
Bleszynski-Jayich et al.,— where aluminium rings were 
used, the high magnetic fields utilized in the experiment 
cause (/) to be negligible, but leave SI unaffectedi^ In- 
deed, the rms fluctuations, given by Eq. ((33|) . agree with 
the measured PCi^ 



VIII. DISCUSSION 

In this work we have studied the disorder-averaged per- 
sistent current of noninteracting electrons. We have ex- 
tended earlier analytical studies, which considered only 
the zero-disorder and the diffusive limitSfi^"— and have 
given an expression, Eq. ([T^. for a general^ ratio of 
L/£, as long as kp£ ^ 1. We define the uncorrelated and 
the correlated-channel regimes in which Eq. (fT3|) can be 
simplified— to the expressions and respectively. 
While previous worksi^"— dealt mostly with ID rings or 
2D cylinders, we have considered here also rings of finite 
narrow width. In particular we have obtained an expres- 
sion for 3D rings. In addition, our expression for the PC 
in a 2D cylinder in the correlated-channel regime in the 
diffusive limit is new. 

The inset of Fig. [2] and Fig. [3] demonstrate that the 
disorder-averaged PC may be a relevant contribution, 
compared with the fluctuation SI, for slightly diffusive 
systems, typically with L/£ < 10. The relation between 

the parameters of a ring that satisfy ( (/m)^ )^^^ > SI, is 
given in Eqs. ([M)) and ([55]) for the uncorrelated and the 
correlated-channel regimes, respectively. We flnd that 
for the parameters of the rings used in Refs. and ITol the 
disorder-averaged PC is relevant compared to SI. 

Interactions, repulsivei^ or attractive can contribute 
to an h/2e flux-periodic disorder-averaged PC. How- 
ever, as long as the sample is not superconducting, 
the PC remains a mesoscopic effect. We have recently 
suggested^^i^ that if the effect of pair-breaking is taken 
into account, attractive interactions can explain the h/2e 
signal measured in ensembles of copper— and goldii 
rings. The contribution of interactions to the PC is 
not sensitive to the exact value of kpL. Therefore, the 
interaction-induced PC may be compared to measure- 
ments using the top equality in Eqs. (j32[) . for any value 
of S{kpL). In contrast, since in reality S{kpL) > tt, 
the interaction-independent contributions of both 51 and 

( (/)^ )^/^ are compared to measurements using the 
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bottom equality in Eqs. ([5^ . Thus, as N increases 
the interaction-dependent contributions to the PC be- 
come dominant over the contributions which do not de- 
pend on electronic interactions. This explains why mea- 
surements on ensembles of 10^ and 10^ rings revealed 
only the h/2e harmonic— It seems that the h/e har- 
monic can be accounted for only by the part of the PC 
that is independent of interactions, which we study here. 
However, since the h/2e periodicity of the interaction- 
dependent part of the PC was obtained from calculations 
of the disorder-averaged PC,— further study is needed 
to assure that the h/e harmonic is not present in the 
interaction-dependent parts of SI. 

Each harmonic has a different temperature depen- 
dence. Higher harmonics decay faster with temperature 
since they necessitate multiple paths around the ring.— 
For this reason we treated the different harmonics sep- 
arately, though our calculations are carried out at zero 
temperature. 

We call attention to the appearance of positive powers 
of the channel number (although the negative power of 
kpL in the correlated-channel regime may partially com- 
pensate that) in the PC magnitude. This implies that 
once multichannel ballistic systems would be manufac- 
tured, relatively large PC's should appear. Both molec- 
ular and clean semiconducting systems come to mind in 
this connection, and perhaps semimetals, such as Bi (see 
first reference of [131). On the other hand, in all regimes, 
the disorder-averaged PC in the diffusive limit is highly 
suppressed by a factor of exp(— L/2£). Again, achieving 
i not too small compared with L, will be helpful. 
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Appendix: An alternative statistical approach for 
the description of the current 

So far we have used the Green function technique for 
our calculations. In this section we develop an alternative 
statistical approach to approximate the current in the 
uncorrelated-channel regime and the zero-disorder limit. 
The following approach leads to the magnitude of the PC, 
which is giveni^ by Eq. (j25p , in a more intuitive way. We 
study here the probabilities that the channels are filled 
with an odd or an even number of electrons, and use the 
results for PCs in canonical ID rings, to obtain the PCs 
of 2D or 3D rings. 



In the regime -1/2 < (/) < 1/2, the PC of a ID ring 
with an odd or with an even number of electrons, see for 
example Ref. 0, is given by 



= [sgn(0) - 2^] 



(A.l) 



(A.2) 



These currents have periodicity of unity in cj). Consider 
a ring of finite width in the grand-canonical ensemble at 
zero temperature. The contribution of the (g, s) channel 
to the PC is obtained by replacing vp hv Eqs. (|A.1|) and 
(|A.2p by an effective Fermi velocity vpiq, s) — Mkp{q, s), 
see Eqs. (jS]) and (fTTj). Here, the exact position where 
the chemical potential crosses the energy levels of each 
channel determines whether the channel is occupied by 
an even or an odd number of electrons, see Fig. 




FIG. 4: The energy levels of a single channel are plotted 
as a function of the flux. The consecutive energy levels for 
a given positive flux and longitudinal indices — n, n, and 
—n — 1, are marked by full circles. The bottom level cor- 
responds to 71 = 0. The random choice of fj, in the interval 
[Eq^s,-n{4'), Eq.s.-n-iifj})] yields an odd number of occupied 
levels when jj, > Eq^s,n{4') and an even number of occupied 
levels when /j. < Eq^s,n {</>)■ The former regime is marked by 
the bold line in the figure. Here, without loss of generality, 
we take n > 0. 

In an ensemble of rings with similar but not identi- 
cal parameters, the energy levels of a given channel are 
shifted (among the rings) due to fluctuations in H and 
W, see Eq. ([4]). Also, the variation of these levels with 
is changing due to fluctuations in L. Therefore, the exact 
position of fi relative to the energy levels of a given chan- 
nel is distributed randomly in the ensemble. When the 
levels with E < Eq^s-n m Fig. |4] are occupied the chan- 
nel consists of an even number of electrons, and when the 
levels with E < Eq^s.n are occupied the channel consists 
of an odd number. Taking Eq^s,7i{<t') — fJ- the probability 
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that a channel consists of an odd number of electrons is 
determined by 



(A.3) 



We assumed here > and n > 0. The difference ap- 
pearing in the nominator is shown in Fig. |4] as a vertical 
line. Inserting the eigenenergies, Eq. in Eq. (jA.3[) 
(considering n 3> 1), yields 



Podd = 1 - 2|0| , P,, 



(A.4) 



These probabilities are independent of the channel index. 

We calculate the average current in an ensemble of sim- 
ilar rings using the currents and the probabilities given 
in Eqs. (|A.ip . \k.2\ . and (|A.4p . and find 



I{q, S) = Podd4dd(g, s) + Pcvcnicvcniq, s) = , 







(A.5) 



For 1 01 1, the probability to have an odd number of 
electrons in a channel is much larger than the probability 
to have an even number, see Eq. (|A.4[) . However, since 
|-?^ovon| > |^dd|, see Egs. ljA.ip and (jA.2p . the average 
PC is zero. This suggests very large fluctuations of the 
current at small flux. The typical magnitude of I{q, s) is 
given by 

1/2 



[l^iq,s)'j = \l fodd/odd(9' «) + ^'ovcn/^vcn('7, «) 

= V2|0|(l-2|0|) . (A.6) 



Here 



Cr 



E 

L q,s 



1 1/2 




(A.9) 



depends on the dimensionality of the ring. The nonana- 
lytic behavior at ^ 1 at zero temperature is due to 
the paramagnetic contributions, since Povcn oc </>, while 
^ovcn oc zbconst at — > 0. Thus, the slope of Eq. (jA.8[) 
at = diverges 4^ 

Equation (jA.8|) reproduces Eq. p5|) obtained for the 
uncorrelated-channel regime in the zero-disorder limit for 
3D rings. For one and two dimensional rings, Eq. (|A.8[) 
reproduces the results of Refs. [l^ andfiol The reason for 
this equivalence is that Eq. ((T8)) . which yields Eq. (|25p . 
is equivalent to Eq. (|A.7p . 

For a finite ensemble of N clean rings, whose typical 
number of channels is A'tot, the probability that all chan- 
nels in all rings will be occupied by an odd number of 
electrons is given for small by 



(Podd) 



NNt, 



ct>NNtat<l 



-> 1 - 2<t>NNtot 



(A.IO) 



We add the assumption that the contributions of differ- 
ent channels to the PC are uncorrelated, which, together 
with Eq. (|A.5p . yields 



I[q,s)I{q',s')^5qq,5ss'P{q,s) 



(A.7) 



Using Eqs. ()A.6j) and (|A.7| we obtain the standard devi- 
ation of the current 



1/2 



(A.8) 



This probability becomes arbitrarily close to unity for 
<j>NNtot ^ 1- Therefore, such a measurement will pro- 
duce the diamagnetic linear response of a clean super- 
conductor (see Sec. |T|. By increasing the flux in a given 
finite ensemble (or by increasing NNtot), even channels 
will appear one by one, each giving a large paramag- 
netic contribution, eventually causing the zero average 
and anomalously large fluctuations of the current. 

Note that an ensemble of ID rings, with equal proba- 
bility for an odd and for an even number of electrons in a 
ring, should exhibit a very large paramagnetic response, 
see Eqs. COj) and (|X2|) . 
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